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A SPECTROSCOPIC SYSTEM EMPLOYING A PLURALITY OF DATA TYPES 
Cross^gferengg to Related Applications 
Thi^PBlication is a continuation-in-part of U.S. patent application Serial No. 
»^^/470^ffledT^^ 22, 1999, which claims priority to and the benefit of U.S. 
provisional patent applica^al No. 60/1 13,761, fded December 23, 1998,^ 
application is related to the U.S^tappHcation entitled, "System for Normalizing Spectra" 
and identified by Attorney Docket NumbWo20, filed on even date herewith, and to the 
U.S. patent application entitled, "Spectral Datamation Of Samples" and identified by 
Attorney Docket Number MDS-02 1 , filed on even da^^h. All of the above applications 
are assigned to the common assignee of this application, and^hereby incorporated by 
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rtnyprnment Rights 

1 This invention was made with government support under Contract No. CA 66481awarded 

bvNauonalCaneerlnstafc.^ 

T?;^ the Invention 
This invention relates generally to spectroscopic systems. More particularly, in one 
embodiment, tine invention relates to a spectroscopic system for characterizing res, specimens 
using a plurality of spectral data types. 

Spectral analysis of biological specimens has been used for disease diagnosis. In general, 
20 spectia a* recorded as va.uesrfanipltode. » measured as a response to an excitation, as 
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a fimction of wavetengti, (or .he inverse of wavetengti, namely frequency). In ft. field of 
sierra, ana ly sis, differen, kinds of information are eonveyed by differen, spectra, types. For 
example, fluorescence speCra are reeorded using a sonree of exci,atio„ illumination ft* ts 
absorbed by a speeimen and «ba, causes the speeimen to emit a fluorescence spec,™ ,ha, 
5 depends in par, on the transfer of energy within and among atoms and/or tnotecutes in the 
specimen. An il.uminarion sonrce for nse in observing and recording flnoreseenee spectra 
g e„era..y operates at a selected monochromatic wavelength, or a narrow range of wave.engtirs. 
Differen. sources of illumination .ha. operate at differen, wavelengths can excite differen, 
, constituents of a specimen. m addition, differen, sources of on Una, opera,e a, differen, 

§ wavetengtins can excite me same constHuen, wi,h differen. efficiencies. Thus, a fluorescence 
t spectr unr can depend on both ,he exemption wavelength tha, is used ,o illuminate a specimen and 
8 „„ the composition of fine specimen itseif. Other effects a.so play a role in determining a 
l fluorescence spectium, for example, ins— effects, effec. relating ,o P o,arization of tire 

S illumination, or thermal effects. 

Some progress has been made in using various optical spectral methods for tura.ysis of 
test specimens, including bio.ogica. specimens. However, the wide varied of physical and 
chemical influences presen, in a tea, speeimen ma, play roles in determining an observed 
spectrnm mate dffficu,, bom tine choice of a suUabte illumination source, and tire interpret 
of the resulting spec*™. Tins is in par, tirte because .here are so many influences on me kind 

effec, relationships among tire multitude of competing influences. 

Cmmarv of lh> Invention 
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The invention provides methods of determining the disease state of a biological specimen 
ohtainedj^r^^^ 

healthy tissue. 

^menfltatisjudgedtobe —ate in creation bv » <~ spec*™ — 
i .odncesafluoreseencespeclramtharisnorcharacrertsticofahcaldrvtiasne. Similar,, 

T.e method comprises recording,*,* flnorescenee and reflectance spectra flom spectflc 
specimen. Themethod then comprises identifying a ptoran* of *P— having a 
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fcamms of the — spectmm residua,. For P«P~ - - 3 ' " 3 

and/or the patient from whom it was obtained. 

^ous.rssue.maap.as fc *-*---^«--*-'««*-- a - 
n ^^iC^eHa.neop.as^CNMn). In anorher em— , medrods of ,ne 
| o^aa^o^—on^the^spec^,™. 
h — — — .ofhefloo.escencespec^aoa.e 

^ — ^««^•--* fc "- ,ta, * fc,- • ,d,,, * 
„»e embodiment, the addiUona. opto, rnformauon is video in—. It. another 
embodiment, the additional optical information is an optical image. 

.es.pecimen.TI.e.s.emcomptisesaaa.collectionmodnletha.collectsafluo.scence 

w nf the nluralitv and observes a reflectance spectrum 
a reflectance spectrum from each member of the plurality, an 

• that is observed to produce a second fluorescence spectrum that is not 
from a test specimen that is observed 10 piu 

a-*' n The svstem further comprises a computation module that 
characteristic of the known condition. The system tunn 
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computes an average reflectance spectrum based npon each member of the phrralitv of firs, 
specimens, and that eonrpu.es a reflectance spectntm residua, by subtracting the average 
reflectance spectrum from the reflect spectrum obtained fiom a test specimen, and an 
analysis module tba, determines the condition of me test specimen based a, leas, in part upon the 

reflectance spectrum residual. 

The foregoing and o,her objecs, aspect feaurres, and advances of ,he invention 
become more apparent from «he following description and from the claims. 
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Brief^gscriEtion " f the drawings 
* „ bj ec,s and features of the invention can be bene, understood with reference ,o .he 
drawings described below, and «be Cairns. The drawings are not necessarUy to scaie, e mP basis 
instead g e„eraU y being placed npon illustrating the princip.es of.be invention. In the drawings, 
l ik . numerals are used to indicate like parts throughout the various views. 

PIG. , shows an exempt spectroscopic system that empioys a plurality of spectra, 
types according .0 an illustrative embodiment of the invention; 

HO. 2 shows an exemplary operational block diagram of the spectroscopic system of 

FIG. 1; 

F,G. 3 is a detailed schematic flow diagram showing exemplary steps of combining a 
„„ce spectrum analysis described more flrlly in F.G. 6 with a reflectance spectrum 
an^sis described more fully in FIG. 7 to perform tissue characterization according ,0 an 

illustrative embodiment of the invention; 

PIG. 4 is a diagram showing an exemplary fluorescence specttam recorded using flte 

spectroscopic system of FIG. 1 ; 

FIG. 5 is a diagram showing an exemplary reflectance spectrum recorded using me 

spectroscopic system of FIG. 1; 

PIG. 6 is a schematic flow diagram depicting an analysis of a fluorescence spectrum in a 

system according to an illustrative embodiment of the invention; and 
20 FIG. 7 is a schematic flow diagram depicting an analysis of a reflectance spectrum in a 

system according to an illustrative embodiment of the invention. 
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rt^tailpH Description 

The invention relates to the use of multiple types of spectral data for determining the 
status of a test specimen In the present invention, a first spectral observation is used to 

5 distinguish between a firs, condition and a plurality of other conditions, and a second observation 
using information obtained from a different type of spectra, dan, is used ,0 distinguish among me 
polity of other conditions. The invention will be described in terms of embodiments tha, relate 
to the use of multiple optical spectra in characterization sysfcms and methods, particulariy in the 
area of medical diagnostics, and especially as such use relates ,o the analysis of spectm obtained 

§ ^hu^ancervicnltissueinthedeWionofcervicnlcancer. However, the invention has 

S applicabilitygenemllyinmeareaofclwac^^ 

UJ plurality of optical spectra types. 

™ KG. 1 depicts an exemplary spectroscopic system 100 employing a plurality of spectral 

| dam types in metitods and systems according to an illuslrative embodiment of tire invention. The 

, 04 accommodates a disposable component 108 which used only once, and discarded after such 
use. The console 102 and the probe 104 are mechanically conneeted by an articulating arm 1 10, 
which can also support the cable 106. The console 102 contains much of the hardware and the 
20 software of the system, and the probe .04 contains the necessary hardware for making suitable 

spectroscopic observations. The details of the system are timber explained in conjunction with 

FIG. 2. 



Atty j^ketNo.: MDS-022 (6219/6) 



8 - 



FIG 2 shows an exempt operational diagmm 200 of a spectroscopic system of 

0 f FlGs . and 2 is sobtanti* .he same as sin,e-beam spec— devices, b», is adapted ,0 
in c,nde te fea t n re sof t he i nve„don.Tneconso,e,02 i nc,udesac Om pn,e r 202 W hichexec^ 

, ^e.ha.connois.heope.donof.hespecnoscopics^.OO. The software inc.udes one 

. ^hme-readabiemedinmcanbe^^ <" 

fcrsoftware.oronecnnsubantu* ftHn^Cc^ins— recordedondevicssocn 
.PROMs, EPROMS oeEEP R OMs,or,he 1 i k e>for software. Tne term machine-readab.e 
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20 ^Mfa^.^-P— ^ rf -"-" ,-,,,,,ieW< " 
chart Of apaften, As described below in more M in an iftostrarive embodimen, of ,he 
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ds eM b,e a user ,0 se,ec, the wavelength rang, for each parfcuiar 
normalization. Other commands enable a user 

^^^^^^^^ , a 
5 ^console.O.alsolncMesanuluavlolet^source.lOsuch.anitrogenlaseror^ 

,. ind „ d e dte ,hecon S o,e,0^oprov M eregu,a,ea.wer f or,eoperat i o„o f ,lo f me 

l oomponen, » — — 

* , ^ooerate with bom the UV light source and the wtute 

W embodiments, a single spectrometer can operate wtt 

? hghrsouree. 

| mdi „3omeemb„dime„ B d ifl e re n t de,ec,orsare»sedmreachlighrsource. ^ 
" ^eonsoMO.alsoine.udeseoupUngoptlcs.O.ocoup.emeUV—.m 

. , irv. rahle 1 06 for transmission to the 

theU V 1 igh.sonree 21 0 t ooneormoreop«c a ,f,bers,nthecable, 

r .he white light illumination from the white light source 212 to one 
probe 104, and for eouphng the whttehght. 

,,Oobservedby.heprobe 104 and carried by one or more oplteal 
on fmm the UV 1 ght source 210 observed oy y 
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observed by the probe 104 and earried by one or more optieal fibers in the cable 106 for 
transmission to the spectrometer and deteetor 218. Tbe eonsole ,02 ine.ndes a foo.swi.eh 224 .0 
enable an opera.or of .he spectroscopic system .00 .o signai when i. is appropria.e to eotnmenee 
a spectral observation by s.epping on .he swi.eh. In .his manner, me opera.or has his or her 
5 hands free to perform other tasks, for example, aligning the probe 104. 

The eonsole 102 inelndes a ealibration port 226 for ealibrating the optieal component of 
.he spee.rome.er system. The operator plaees .he probe .04 in regislry win. .he calibration port 
226 and issues a command ma, sterts the calibration operation. ,n the calibration operation, a 
, , calibrated ligh, sonrce provides il.nmination of known in.enstiy as a taction of wave,eng,h as a 
1 calibration signal. The pnah, 104 de.ec, .he calibration signal. The probe .04 .ransmits me 
S de.ec.ed signal through me optica, fiber in the cable ,06, through the coupling optics 222 ,0 me 
U spec,rometeranddetec.o,218. A test spectral result is obtained. A calibration of the spectral 
l system is computed as .he ratio of me amplimde of the known illumination a, a particular 
1 wavelengm divided by the test spectral resul, at the same wavelengm. 

§ The probe 104 includes probe optics 230 for illuminating a specimen «o be analyzed with 

™ UV and white light from me UV source 2.0 and me white light source 2,2, and for collecting me 
fluorescent and backscatter (or reflectance) illumination from the specimen ma, is being 
a^lyzed. The probe includes a scanner assembly 232 that provides illumination from me UV 
source 2,0 in a raster pa«em over a target area of .he specimen of cervica, tissue ,0 be analyzed. 
20 The probe includes a video camera 234 for observing and recording visual images of .he 

determine where on me surface of the specimen to be analyzed the probe ,04 is pointing. The 
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probe 104 also includes a white light illuminator 238 to assist the operator in visualizing the 
specimen to be analyzed. Once the operator aligns the spectroscopic system rtnd depresses the 
footswitch 224, the computer 202 controls the actions of the Ugh, sources 210, 2.2, the coupling 
optics 220, .he transmission of light signals and electrical signals through the cable .06, the 
5 operation of the probe optics 230 and the scanner assembly 232, the retreiva. of observed spectra 
via the cable 106, the coupling of the observed spectra via the coupling optics 222 into the 
spectrometer and detector 218, the operation of the spectrometer and detector 218, and the 
subsequent signal procesing and analysis of the recorded spectra. 
5 HO. 3 is a detailed schematic flow diagram 300 showing exemplary steps of combining 

I the fluorescence spectrum analysis described more Ml, in HO. 6 below with the reflectance 
« speomrm analysis described more fully in FIO. 7 below ,0 perform tissue characterization 
™ according to an illustrative embodiment of the invention. Step 3.0 indicates that the results 
| presented in FIG. 6 for analysis of fluorescence spectra torn a test specimen of unknown 
1 condition or unknown state of health are available. At step 320, the computer 202 determines 
1 whether the test specimen can be classified as "normal," or "metaplasia," or can no. be classified 
by fluorescence spectroscopy alone. This process is described in delail a. s.ep 665 of FIG. 6 
below. As indica.ed in s.ep 325, a decision is taken as <o whether the tea. specimen has a 
definitive state of heahh, for example «ha, .he specimen is "norma,." If .he .est specimen can be 
classified, for example as normal, the method ends at step 330. 
20 In the event tha, a definitive condition or state of heal* cannot be ascribed to a test 

specimen, the computer 202 further analyses information available from a reflectance spectrum 
or from a plurality of reflectance spectra taken from the test specimen. A. step 335, the computer 
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202 pmvides reflec,ance speCra processed according .0 me systems and methods described in 
connection with FIG. 7 below. 

If the specimen eanno, be eiassifred, a mean nortnabzation step is performed by computer 
202 as indicate* a, step 340. Tbis mean „orma,iza,ion step is described in detai. a, step 755 of 

u speCra that indicate tba, those locations can be Cassified as representing nomra, smramoos 

1 4a, are no, capabie of being classified as "norma," or »me»p,asia» so,e, y on me basis of 

i y 

i« fluorescence spectra. 

S AS indica,ed in step 350, ,he compn,cr 202 ean carry on, an analysis nsing a metric as 

m e tri e i nN.dimensi„na 1 space. .» one embodiment, the test reflectance spec*a me trttncafcd ,o 
the wavCength regions 39, nm ,o 4S4 nm, and 532 nm ,o 625 nm. In one embodiment, the 
creations OH , and CIN flfll are the Cassations ma, are possib.e for a test spec™, tba, 

202 provided ,ba, ,he va,ne of me metric does no, exceed a pre-de.erntined maximnm va,ne. 
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■ At step 360, the computer 202 presents 'the results of the classification of the test 
specimen, as a condition or state of heaiO, corresponding to one of norma,, metaplasia, ON I and 

cm n/ra. 

FIG. 4 is a diagram 400 showing an exemplary fluorescence spectrum recorded using the 
5 spectroscopic system of FIG. 1. In FIG. 4, a curve 4.0 having an amplitude defined in terms of 
number of counts, for example a signal strength expressed as numbers of photons, as indicated 
along the vertical axis 420. An amplitude is plotted against a wave.ength of light expressed in 
nm, as indicated along the horizontal axis 430. Such spectra can be recorded using the methods 
3 . , and systems described in Figs. 1-3 with an illumination source that provides optical excitation 
M comprising a substantia,* monochromatic source in the n„rav,o,e, portion of the optica, 

„ . , m ^pnjpted in FIG 4 was recorded from human 

spectrum. The exemplary fluorescence spectrum depicted m riu. * 

UJ 

m cervical tissue using ultraviolet illumination. 

?, FIG. 5 is a diagram 500 showing an exemplary reflectance spectrum recorded using the 

1 specuoscopiesystemofFIG.1. ,nF,G. 5, a curve 5.0 having an amplitude defined in terms of 
« number of counts, for example a sigaa. stiength expressed as numbers of photons, as indicated 
a,ong the vertiea, axis 520. A* amphmde is p.otted against a wave.ength of Ugh. expressed in 

as indicated a,o„g flte horizonta, axis 530. Such specUa can be recorded using the methods 
and systems described in Figs. ,-3 with an inumination source that provides optica, excitation 
comprisingabroadbandson^einthevisib.eportionof.heop.ica.spectnnn. Theexemp.ao- 

, • x j • ct/- * ,1,3c rprnrded from human cervical tissue using 
20 reflectance spectrum depicted in FIG. 5 was recorded irom nu 

broadband visible light illumination. 
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FIG. 6 is a schematic flow diagram 600 depicting an analysis of a fluorescence spectrum 
in a system according ,0 an illustrative embodiment of the invention. Fluorescence spectra sod- 
as the exemplary spec.™, of FIG. 4 are suitable for the analysis described hereafter. The 
spectral analysis involves the measurement of spectra from a res. specimen whose condition is to 
be determined, and employs a plurality of reference spectra taken from one or more specimens 
having known conditions. Those of ordinary skill in the spectroscopic arts will underhand that 
the reference spectra can be obtained befote, after, or a, the same operation session as the test 
spectrum or spectra are obtained. For purposes of exposition, .he treatment of .he reference 
spectra will be described first. The treatment of a test spectrum will .hen be described. 

In the present application, the term "characteristic N-dimensional value" should be 
understood to comprehend whichever of a point, a volume, a surface, a probability distribution or 
the like is used in the subsequent analysis. In mathematical tetms, the "characteristic N- 
dimensional value" can be understood ,„ be an ordered N-tuple of values, each one of which 
values can be expressed in its own dimensional units. It is well known in the spectroscopic arts 
to record one or more reference spectra from specimens that are known to have specific 
conditions. For example, it is known to record reference spectra from specimens having known 
states of health or known disease conditions. Reference spectra can be manipulated using the 
same pre-processing methods tha, are applied to spectra recorded from tea, specimens, so as to be 
able to compare spectra having substantially similar processing histories. 

The compote, 202 can store in a machine-readable memory for later use .he information 
, comprising .he various characteristic N-dimensio„al values that are computed, as well as the 
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corresponding known condition. The compote; 202 can also similarly store in a machine- 
readable memory the various reference reflectance spectra data to were obtained. 

Turning to the discussion of the recording and analysis of a test spectrum, as indicated in 
FIG. 6 a, step 610, an operator collects and records a test fluorescence speetntm or a plurality of 
5 test fluorescence spectra from a test specimen of unknown state of heahh or unknown condition. 

S,e P 615 indicates that the computer 202 can smooth the raw data. Various mathematiea! 
smoothing techniques are known in the art, such as computing moving averages, or applying 
filters to remove noise, for example using well-known convolution methods. 

Step 620 indicates that the computer 202 con subtract a background reading, for example 
g a reading taken using the systems and meurods of FIOs. .-3 with a standard specimen, or taken 
2 without a tes, specimen, from the reading observed from a test specimen or a reference specimen. 
S ,„ step 625, fine computer 202 can average a plurality of spectral data taken from the same 
'I specimen, for example, the computer 202 can average N repeated spectral measurements taken 
I from the same location of a single specimen, where N is an integer greater than 1 . 
1 In step 630, the computer 202 can correct the wavelength assigned to a particular 

m amplitudeofaspecnw.orcaneo^.aplura^ofsuchwave.engflts. The computer 202 can 
perform such a correction by using a pre-recorded calibration file recorded from an optica, source 
having fines a. known wave.engflts, for example a spectrum recorded from a lamp having the 
chamcteristic emission fines of mercury, or another we.l knotvn emission source. The computer 
20 202 uses the known wavelengfits of the characteristic emission lines ,„ detenuine a relationship 
between wavelength and a data channel number assigned by the computer 202 and an A/D 
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converter to the characteristic spectra. The computer 202 cao determine a wavelength for each 
A/D channel and thus knows the correct wavelength to assign to each amplitude. 

As indicated in step 635, the computer 202 can correct the spectral data for effects 
ascribable to open air, such as absorption due to atmospheric constituents, for example gases 
such as water vapor, carbon dioxide, and others mat have strong absorption a, characteristic 
wavelengths. 

Step 640 indicates mat the computer 202 can correct for instrument-induced features by 
dividing an observed spectrum by an instrument function. The instrument function is obtained 
by dividing a known, accurate spectrum, such as a NIST-traceable tungsten lamp spectrum, by 
the observed spectrum recorded for illumination from such a lamp passing through the 
instrument of FIGs. 1 and 2. Dividing a firs, spectrum by a second spectrum involves dividing 
the respective amplitudes of the firs, spectrum by the corresponding amplitudes for the second 
spectrum and recording theresul. of each such division in a file as a function of the respective 
wavelengths, such as the instrument function. 

Step 645 indicates that the computer 202 can truncate a spectrum to limit tire data saved 
,„ a file as amplitudes corresponding to the reduced set of wavelengths, rather than all the 
wavelengtits that the spectral instrument of FIGs. 1-2 is capable of recording. The truncation is 
performed to avoid saving data that has no significance, such as data corresponding to amplitudes 
having zero intensity for wavelengtits beyond some predetermined value. In some embodiments, 
) the operator can use the keyboard 206, or another input mechanism, to indicate one or both 
wavelength limits of the truncated wavelength range to be retained. 
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As indicated at step 650, the computer 202 can reduce the number of amplitudes, each 
amplitude corresponding to a selected wavelength, that are used to characterize a spectrum. In 
some embodiments, the number of amplitudes used to characterize a spectrum is fifty (50). In 
principle, as few as three amplitudes at three wavelengths may be sufficient to characterize a 
fluorescence spectrum as representing a test specimen comprising healthy normal squamous cells 
as distinguished from other tissue types having other states of health. 

As indicated at step 655, the computer 202 can normalize the reduced spectral data using 
a system and method called normalization using non-uniform segmentation. Normalization 
using non-uniform segment normalization is described in detail in the co-pending patent 
application entitled "System for Normalizing Spectra," which application is commonly assigned 
to the assignee of the present application, and which application is incorporated herein by 



!*j reference in its entirety 
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In one embodiment, where the metric to be used in classifying a test specimen is a 
1 Mahalanaobis distance, the computer 202 performs a mathematical procedure intended to 
1 guarantee that a matrix used in former analytical steps is capable of being inverted. Matrix 
inversion is a mathematical process well known in the matrix mathematical arts. In one 
embodiment, the matrix that is computed is called a Friedman matrix. The data used in the 
calculation of the Friedman matrix is obtained from reference data recorded and stored in 
machine-readable format. The step of computing the appropriate Friedman matrix can be 
20 performed at any time after the reference spectral data is available. 

Computation of the Mahalanobis distance requires the inversion of a weighting matrix. 
One approach (i.e., Linear Discriminant Analysis) utilizes the inverse of the pooled within- 
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groups covariance while another method (Quadratic Discriminant Analysis) employs the inverse 
of each within-group covariance. When a "large" number of wavelengths are used (e.g., for large 
p), either of the aforementioned matrices can be singular and hence cannot be inverted. 
Therefore, an alternate weighting matrix called the Friedman matrix is used. 

7 The Friedman matrix is the weighted linear combination of a within-group covariance 
matrix Cj and the pooled within-group covariance C that is a function of the Friedman 
parameters y and X. These parameters are selected from the unit interval [0,1]. It is important to 
note that the Friedman parameters are not physical quantities whose values carry classification 
rformation. Rather, they are used to insure that the weighting (Friedman) matrix is non- 
ngular. Equation 1 details the p-by-p Friedman matrix, which is denoted by XI 

Q.^-O-rt^-^-Ci^-^^-^^^^"* 7 - Eqn ' 1 
N note that *A) is the trace of the matrix A, namely the sum of the diagonal elements. Also, 
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/ Yn is the p-dimensional identity matrix. 
pxp 
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Let s(A p ) be ap-dimensional spectrum. The Mahalanobis distance d from s(A p ) to the 
group mean u; with respect to the associated Friedman matrix Ufi, X) and Friedman parameters 
y and X is given in equation 4 below. 



5 4(A p ),M l )=^iA r )-M J )'a i (r,A)-'^)-M i r\ E «"' 4 

I 

Where other metrics are used to classify a test spectrum and a test specimen, the step 670 can be 

replaced by the appropriate computation. 
| As indicated at step 660, the computer 202 computes a metric from each test spectrum (or 

| from an average value representing a plurality of spectra, such as performed at step 625) to each 
S characteristic N-dimensional value for the reference spectra recorded. In one embodiment, the 
™ metric can be a Mahalanobis distance, as indicated above. In other embodiments, the metric can 
1 be the square root of the sum of the squares of the differences in coordinates in the N- 
I dimensional space, or the metric can be the Bhattacharyya distance. 

H At step 665, the computer 202 examines the results obtained in step 660, and determines 

a classification for the test specimen, based on the metrics computed with respect to the test 
fluorescence spectrum and a classification rule or relation. In one embodiment, the classification 
relation is that the condition of the test specimen is assigned as the condition corresponding to 
the reference spectrum constellation or set having the Mahalanobis shortest distance, provided 

20 that the shortest Mahalanobis distance is less than a predetermined minimum distance. In one 
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embodiment if no M*a.a„obis distance is less" than a predetermined minimum distance, me test 
spectrum is discarded as being incapable of being classified, or indeterminate. 

,„ one embodiment, the tes, fluorescence spectrum can be compared to multiple sets of 
reference spectra in a sin gl e comparison, .n one embodiment, the computations can be repeated 
for different sets of reference conditions. For example, a res. fluorescence spectrum from a test 
specimen can be compared in a flrs. computation to reference fluorescence spectral data for 
norma, squamous tissue and CIN Mil tissue, compared in a second computation to reference 
fluorescence spectra, data for nornra. squamous tissue and CDS . tissue, and compared in a third 
computation to reference fluorescence spectra, data for norma, squamous ttssue and metaplasia 
.issue. If the tes, fluorescence spectrum is classified in each of the three comparisons as being 
more closely related to the reference spectra for notma. squamous tissue, the test specimen can 
be classified by the computer 202 as having normal health, and the classification process is 
complete. However, if the test fluorescence spectrum is classified in either of the comparisons 
involving CIN I or CIN Mil as being less closely related to the reference spectta for normal 
squamous tissue man one of .he other conditions or states of health (namely ON ,, or ON Mil), 
me computer 202 can report that the state of heal* of the tes, specimen is no, clear.y ,ha, of 
normal squamous tissue, and ma, further analysis is in order. In one embodimen,, tire further 
aralysis involves the recording and examination of reflectance spectta. 

FIG. 7 is a schematic flow diagram 700 depicting an analysis of a reflectance specttum in 
0 a system according ,0 an illustrative embodiment of the invention. RefleOance spectta such as 
me exemplary specttum of F.C-. 5 are suitable for the analysis described hereafter. The spectta. 
analysis involves ,he measurement of spectta from a ,es, specimen whose condition is ,o be 
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spectrum or spectra are obtained. 

l nG 7 ^ W -*----— — — — — °" Pta *° f 

, ft™ a, est specimen of unknown state of health or unknown condmon. 
sO test reflectance spectra from a test specimen 

g Step715ind ,ca t es,har ft ecompo,er 2 02cansmoo,h a ,erawaa.. Various mathematrca! 

A f v.mnleusine well-known convolution methods. 

5 filters to remove noise, for example using wen 

« ♦ onor^n subtract a background reading, for example 

£ step 720 indicates that the computer 202 can subtract an g 

, mpthods 0 f FIGs 1-3 with a standard specimen, or taken 
a reading taken using the systems and methods ol Ml* 

• ^mthereadingobservedfromatestspecimenorareferencespecimen. 
without a test specimen, from the reading ods 

^.epm.ecomputer^canavera.eapiuraUtyofspec^o.tarakenfrommesame 

„ step 730, the computer 202 can correct the wave,en g th assigned to a particuia, 
^.udeofaspectrum.orcancorrectap.umh.ofsuchw— , The computer 202 can 
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perform such a correction by using a pre-recorded calibration file recorded from an optical source 
having lines a. known wave.eng.hs, for example a spectrum recorded from a lamp having the 
characteristic emission lines of mercury, or another wel. known emission source. The computer 
202 uses the known wavelengths of the characteristic emission lines ,o determine a relationship 

5 between wavelength and a data channel number assigned by the computer 202 and an A/D 
converter to the characteristic spectra. The computer 202 can determine a wavelength for each 
A/D channel and thus knows the correct wavelength to assign to each amplitude. 

As indicated in step 735, the computer 202 can correct the spectral data for effects 

* ascribable ,0 open air, such as absorption due ,0 atmospheric constituents, for example gases 

S such as water vapor, carbon dioxide, and others that have strong absorption a, characteristic 

SI 

^ wavelengths. 

8 step 740 indicates una. the computer 202 can correct for instrument-induced features by 

l dividtng an observed spectrum by an mstrume„, tenon. The instrument function is obtained 
| by dividing a known, accurate spectrum, such as a NIST-traceab.e tungsten .amp spectrum, by 
,| the observed spectrum recorded for illumination from such a lamp passing through the 

instrument ofFICJs. 1 and 2. Dividing a firs, spectrum by a second spectrum involve, dividing 
the respective amp.itudes of the first spectrum by the corresponding amplitudes for fire second 
speetrnm and recording the result of each such diviaton in a file as a function of the respective 
wavelengths, such as the instrument function. 

Step 745 indicates that the computer 202 can truncate a spectrum to limit the data saved 
file as amplitudes corresponding to fire reduced set of wavelengrhs, rather than all the 
wavelengths that the spectra, instrument of FIG, 1-2 is capable of recording. The tnrncation is 
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performed to avoid saving data that has no significance, such as data corresponding to amplitudes 
having zero intensity for wavelengths beyond some predetermined value. In some embodiments, 
the operator can use the keyboard 206, or another input mechanism, to indicate one or both 
wavelength limits of the truncated wavelength range to be retained. 

As indicated at step 750, the computer 202 can reduce the number of amplitudes, each 
amplitude corresponding to a selected wavelength, that are used to characterize a spectrum. In 
some embodiments, the number of amplitudes used to characterize a spectrum is fifty (50). In 
principle, as few as three amplitudes at three wavelengths may be sufficient to characterize a 
reflectance spectrum as representing a test specimen comprising healthy normal squamous cells 
S as distinguished from other tissue types having other states of health. 

As indicated at step 755, the computer 202 can normalize the reduced spectral data using 
f a system and method called mean normalization. Mean normalization involves determining a 
1 mean value at each wavelength of interest for reference reflectance spectra and subtracting the 
I mean values so determined at a particular wavelength of interest from the amplitude of the test 
ffi spectrum (or spectra) at the same wavelength. 

In one embodiment, at step 760, where the metric to be used in classifying a test specimen 
is a Mahalanaobis distance, the computer 202 performs a mathematical procedure intended to 
guarantee that a matrix used in further analytical steps is capable of being inverted. Matrix 
inversion is a mathematical process well known in the matrix mathematical arts. In one 
20 embodiment, the matrix that is computed is called a Friedman matrix. The data used in the 

calculation of the Friedman matrix is obtained from the stored data recorded in machine-readable 
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format. The step of computing the appropriate Friedman matrix can be performed at any time 
after the reference spectral data is available. 

Computation of the Mahalanobis distance requires the inversion of a weighting matrix. 
One approach (i.e., Linear Discriminant Analysis) utilizes the inverse of the pooled within- 
groups covariance while another method (Quadratic Discriminant Analysis) employs the inverse 
of each within-group covariance. When a "large" number of wavelengths are used (e.g., for large 
P ), either of the aforementioned matrices can be singular and hence cannot be inverted. 
Therefore, an alternate weighting matrix called the Friedman matrix is used. 

The Friedman matrix is the weighted linear combination of a within-group covariance 
matrix Cj and the pooled within-group covariance C that is a function of the Friedman 

§ parameters^. Tta-p-^-e^^^^^^ * " * 
note that the Friedman parameters are not physical auantities ,hose values carry classification 
information. Rather, they are used to insure that the weighting (Friedman) matrix is non- 
1 singular. Equation 5 details thep-by-p Friedman matrix, which is denoted by'Q. 

Q /(ri A) = (l-rt-kl^)-C j ^-cl^'»-[0-^-C;^-c]-/„ E,n.5 

N note to rr(A) is .he trace of the matrix A, namely the sum of the diagonal elements. Also, 
I pxp is the p-dimensional identity matrix. 
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Eqn. 7 
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Let S (h p ) be ap-dimensional s P ec,rum ' The Uahalamhis d>s ' mC " d fr ° m S< V '° "* 
group mean p, with respec, to the associated Friedman matrix fi/y, X) and Friedman parameters 

Y and X is given in equation 8 below. 



Eqn. 8 



Where other metrics are 



used to classify a test spectrum and a test specimen, the step 770 can be 



replaced by the appropriate computation. 

As indicated at step 760, the computer 202 computes a metric from each test spectrum (or 
from an average value representing a plurality of spectra, such as performed at step 725) to each 
characteristic N-dime„sional value for .he reference spectra recorded. In one embodiment, the 
metric can be a Mahalanobis distance, as indicated above. In other embodiments, the metric can 
be the square root of the sum of the squares of the differences in coordinates in the N- 
dimensional space, or the metric can be the Bhattacharyya distance. 

At step 765, the computer 202 examines the results obtained in step 760, and detemtines 
a classification for the test specimen, based on the metrics computed with respec, to the test 
reflectance spectrum and a classification rule or relation. In one embodiment, the classification 
relation is that the condition of the test specimen is assigned as the condition corresponding to 
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. U tinn or set having the Mahalanobis shortest distance, provided 
the reference spectrum constellate or set having 

. then a nre-determined minimum distance. In one 
that the shortest Mahalanobis distance is less than a pre determ 

w a- . nne is less than a predetermined minimum distance, the test 
embodiment, if no Mahalanobis distance is less than p 

fordif ferentsetsofreferencecondition S . F or example, a test re fl ectance s P ect,m from a test 
sp „be^^ 

1 dby 

. • omnal health and the classification process is complete. However, 
the computer 202 as having normal health, ana 

, ™ l ™ ON Mm, the computer 202 can report 
fl.eothercondWor.orsta.esofhealthlname.yCINl.orCINlWl ), 

20 examination of reflectance spectra. 
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Equivalents 

without departing from the spirit and scope of the invention 



form and detail may be made therein 
5 as defined by the appended claims. 

What is claimed is: 
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